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<^ . 1. Introduction 

. The purpose of the current heavy ion programs at CERN (Switzerland) and Brookhaven 
National Laboratory (USA) is to probe strongly interacting matter under extreme con- 
ditions, i.e. at high densities and temperatures. The central subject of these studies is 
the transition from the quark-gluon plasma to hadronic matter. In the early phases of 
ultra-relativistic heavy ion collisions, when a hot and dense region is formed in the center 
of the reaction, there is copious production of up, down, and strange quarks. Transverse 
expansion is driven by the multiple scattering among the incoming and produced parti- 
O \ cles. As the medium expands and cools, the quarks/gluons combine to form the hadrons 
^ that are eventually observed. 

In June 2000, the long awaited gold collisions were delivered by the RHIC collider 
at -Js NN = 130 GeV. A great amount of new data, analyzed from that brief period of 



running, was presented at this conference |j]-f|]. From SPS, results on the systematic of 
the transverse momentum distributions of <f> and J /if) ||[7| were shown. New results 
from 40 and 80 A- GeV (a/s nn ~ 9 and 12 GeV) collisions were also shown. 
^ ■ In this paper, we will summarize the recent experimental results on transverse momen- 

tum distributions and particle ratios. The related physics issues are chemical equilibrium 
and collective expansion. While the former has to do with the inelastic collisions, the 
later is dominated by the elastic cross section. We will focus at the relatively low trans- 
verse momentum region, p t < 2 GeV/c, where the bulk production occurs. For results on 
high momentum transfer, global measurements, and correlations readers are referred to 
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2. Kinetic freeze-out: transverse momentum distributions 

At this conference, identified particle transverse momentum distributions at mid-rapidity 
were shown by the BRAHMS, PHENIX, and STAR experiments 0,|,[§|]]||g. While 
PHENIX measured tt, K, and p up to 2.5 GeV/c in p t , the STAR TPC provided clean 
spectra of identified particles up to about 1 GeV/c in p t within a somewhat wider rapidity 
window. The preliminary transverse momentum distributions for negative pions, kaons, 
and anti-protons from the STAR Collaboration are shown in Fig. 1. The results are from 
the 6% most central Au+Au collisions. Hydrodynamic motivated fits are also shown in the 



figure as dashed lines [15]. The extracted common freeze-out temperature and collective 
velocity parameters are about 100 MeV and 0.6c, respectively. In Fig. 2, the preliminary 
transverse momentum distribution of anti-protons from the PHENIX collaboration |l3j 
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is compared to the STAR results [14]. It is important to note that the transverse mo- 
mentum distributions are consistent within the overlap region. Below m t — mass = 0.7 
GeV/c 2 , both sets of data are consistent with the hydrodynamic fit. The lower part of 



the distribution is from the NA44 Collaboration at SPS 



17.2 GeV) 
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Figure 1 The STAR preliminary results 
of the transverse mass distributions for 
negative pions, kaons, and anti-protons 
from Au+Au central collisions. Dashed 
lines represent the model fit results. 
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Figure 2 PrelintytoaVffi^k^Pofdfr trans- 
verse mass distributions from STAR and 
PHENIX (circles) are compared. Squares 
represent the results from central Pb+Pb 
collisions. Open symbols represent the 
model fit results. 



The measured transverse momentum distributions also have been fitted by the ex- 
ponential function f = A ■ exp(—m t /T), where T is the slope parameter and A is the 
normalization constant. The magnitude of the slope parameter provides information on 
temperature (random motion in local rest frame) and collective transverse flow. Fig. 3(a) 
shows the measured particle slope parameters from Pb+Pb central collisions at SPS (a/s nn 



8.8 GeV and 



17.2 GeV) energies []8|,p|, [T6| , [TT|1 . The new results are: 



Particle distributions from 40 A- GeV (y/s NN = 8.8 GeV) collisions were reported 
by NA45 and NA49 HJJ. The slope parameters are similar to the ones observed at 
158 A- GeV, 



i.e.. 



they follow the established systematic trend [18 



For the first time, the systematic of the transverse momentum distributions for 
charm particles (J /if?) was reported by the NA50 experiment |||7|]. It is interesting 
to note that the slope parameter of J/ip is similar to those for and Q; 

For central collisions the <fi slope parameter of NA49 is about 300 MeV |J whereas 
that of the NA50 is about 240 MeV. NA49 and NA50 reconstructed <p mesons 
via K + K~ and channels, respectively. As discussed in |l9| , [20|[ , part of the 

difference may be caused by the final state interaction of the decay kaons. In fact, 
if one studies the centrality dependence of the <fi slope parameter, one finds that at 
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peripheral collisions, both experimental results agree with each other and the value 
is close to that from p + p collisions . 
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Figure 3 Slope parameters as a function of particle mass for (a) Pb+Pb central colli- 
sions at SPS (^/s NN = 17.2 GeV) and Au+Au central collisions at (b) RHIC (y/s NN = 
130 GeV). Weak and strong interacting limits are indicated by I and II, respectively. 

Preliminary results of the slope parameters from Au+Au central collisions at RHIC 
(\/ s nn = 130 GeV) [|],[14] are shown in Fig. 3(b). It is obvious that the mass dependence 
of the slope parameter is stronger than that from collisions at SPS energies. 

As one can see from Fig. 3(a), the slope parameters appear to fall into two groups: 
group (I) is flat as a function of particle mass, whereas the slopes in group (II) increase 
strongly with particle mass. At the RHIC energy, the slope parameter systematic of 
7r, K, and p shows an even stronger dependence on the particle mass. The strong energy 
dependence of the slope parameter might be the result of the larger pressure gradient at 
the RHIC energy. With a set of reasonable initial/freeze-out conditions and equation of 
state, the stronger transverse expansion at the RHIC energy was, in deed, predicted by 
hydrodynamic calculations [21-23]- In addition, the results are consistent with the large 
value of the event anisotropy parameter v 2 at the RHIC energy P3||?5| . 

Within the hadronic gas, the interaction cross section for particles like <p, Q, and J/i[) are 
smaller than that of ir, K, and p [ 26| . Therefore the interactions between them and the rest 



of the system are weak, leading to the flat band behavior in Fig. 3(a). On the other hand, 
the slope parameter of these weakly interacting particles may reflect some characteristics 
of the system at hadronization. Then it should be sensitive to the strength of the color 
field " 
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Under this assumption, the fact that the weak interacting particles show a 
flat slope parameter as a function of their mass would indicate that the flow develops at 
a later stage of the collision. Should the collective flow develop at the partonic level, one 
would expect a mass dependence of the slope parameters for all particles |p0 |. 

Based on the idea of Color Glass Condensate, McLerran and Schaffner-Bielich pre- 



dicted |H| that the measured mean p t is controlled by the intrinsic p t broadening rather 
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than transverse flow. It will be interesting to see the slope parameters of 4>,Q, and J/ip 
from collisions at RHIC energies. As indicated by the horizontal bar in Fig. 3(b), one 
might expect both the absolute magnitude and the slope as a function of particle mass to 
be different from that observed at SPS energies. 

Figure 4 shows the bombarding energy dependence of the freeze-out temperature Tf 
and the average collective velocity (f3 t ). It is interesting to observe that both Tf a and 
(Pt) saturate at a beam energy of about 10 A-GeV. The saturation temperature is about 
100 - 120 MeV, very close to the mass of the lightest meson [32]. The steep rise of the 
Tf and ((3 t ) up to about 5 A- GeV incident energy indicates that at low energy collisions 
the thermal energy essentially goes into kinetic degrees of freedom. The saturation at 
~10 A-GeV shows that particle generation becomes important. As proposed in [ |32"| , |33" 



for a pure hadronic scenario there may be a limiting temperature T c fa 140 MeV in high- 
energy collisions, although the underlying physics for both, the transition from partonic to 
hadronic degrees of freedom and the transition from interacting hadrons to free-streaming 
is not clear at the moment. By coupling the limiting temperature idea to a hydrodynamic 
model calculation, Stocker et al. successfully predicted |34] the energy dependence of 
the freeze-out temperature. It is worth noting that measurements of the pion phase 
density p5|-p^1 at freeze-out also show saturation at E beam ~10 A-GeV. 
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Figure 4 Systematics of kinetic freeze-out temperature parameter Tf and average collective 
transverse flow velocity (fit) as a function of beam energy. At ^s NN ~ 5 GeV, both values 
of temperature and velocity parameters seem to saturate. However, the velocity parameter 
extracted from central collisions at the RHIC energy is higher than values from collisions at 
lower beam energies. 

At the RHIC energy, the collective velocity parameter seems to be larger than that from 
collisions at AGS/SPS energies. This can already be seen in Fig. 3 where the increase from 
pion to proton at y/s NN = 130 GeV (Fig. 3(b)) is much faster than at y/s NN = 17 GeV 
(Fig. 3(a)). Is this the manifestation of van Hove's [0 picture? On the other hand, 



compared to results from lower energy collisions, the temperature parameters seem to be 



lower. Is this the consequence of hydrodynamic expansion from a higher initial density 
fireball? An energy scan between y/s NN = 20 - 130 GeV will be extremely important in 
order to study this evolution in more detail. 
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3. Mid-rapidity p/p ratios and net-proton distributions 

All RHIC experiments have consistently measured the ratio of p/p at mid-rapidity 
ED|-H. The dependence of the p/p ratios as a function of the center of mass energy is 
shown in Fig. 5. Only central collision results were used for heavy ion collisions. Clearly, 
although the ratio is not unity, a dramatic increase in the ratio is observed from SPS to 
RHIC, meaning that the system created at a/s nn = 130 GeV is close to net-baryon free. 
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Figure 5 Mid-rapidity p/p ratios measured 
in central heavy-ion collisions (filled symbols) 
andp + p collisions (open symbols). The left 
end of the abscissa is the p-p pair production 



threshold in p + p (. 



-3.75 GeV). 
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Figure 6 Net-protons yields as a function 
of charge particle multiplicity. Open and 
filled symbols represent the yields within 0.3 
< Pt < 1-0 GeV/c and the values extracted 
from the full transverse momentum range, 
respectively. No hyperon decay corrections 
have been applied to the yields. 



One of the unsolved problems in high energy nuclear collisions is the baryon transfer 



that occurs at the early stage of the collision [ |5| . The later dynamic evolution of the sys- 
tem is largely determined at this moment since the available energy is fixed at this stage 
of the collision. On the theoretical side, a novel mechanism, the non-perturbative gluon 
junction, was proposed to address this problem [|6|,|7|. From the STAR measured p/p 
ratios and the STAR preliminary results of anti-proton yields [|j], we have extracted 
the net-proton yields at mid-rapidity. The values are shown in Fig. 6 as a function of 
collision centrality. Open symbols represent the measured yields within \y\ < 0.1 and 0.3 
< Pt < 1-0 GeV/c. The fill symbols show the values extracted from the full transverse mo- 
mentum range. Estimated systematic uncertainties are shown as caps. No hyperon decay 
corrections have been applied as that requires the knowledge of the yields of hyperons. 
The hatched bar represents the BRAHMS preliminary results [l|]. 

One can see that the number of net protons increase as the number of the negatively 
charged particles increases. This implies that more and more protons (baryons) are trans- 
fered from the incoming nuclei to mid-rapidity as the overlap of the two nuclei increases. 



More data are needed to study the role of the junction mechanism [|46| , f47| in heavy ion 
collisions. In this respect, the rapidity distributions of the net-protons as a function of 
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collision centrality will be of crucial importance. 



4. Chemical freeze-out: particle ratios 

Recently, much theoretical effort has been devoted to the analysis of particle production 
within the framework of statistical models []4g| , |49|] . These approaches are applied to the 
results of both elementary collisions (e + e~ , p + p) and heavy ion collisions (Au + Au and 
Pb+Pb) |j48"| . Many features of the data imply that a large degree of chemical equilibration 
may be reached both at AGS and SPS energies. The three most important results are: (i) 
at high energy collisions the chemical freeze-out (inelastic collisions cease) occurs at about 
160-180 MeV and it is 'universal' to both elementary and heavy ion collisions; (ii) the 
kinetic freeze-out (elastic scatterings cease) occurs at a lower temperature ~ 120 — 140 
MeV; (iii) the compilation of freeze-out parameters ]50| in heavy ion collisions in the 
energy range from 1 - 200 A-GeV shows that a constant energy per particle (E) / (N) ~ 1 
GeV can reproduce the behavior in the temperature-potential (T c h — fi B ) plane fl50| . 
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Figure 7 STAR preliminary particle 
ratios vs. results of thermal model fits. 
The thermal parameters for tempera- 
ture and chemical potential are T c h = 
190 ± 20 Me V and fi B =45± 15 MeV, 
respectively. 
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Figure 8 Phase plot T c \ t vs. fi B . Dashed-line 
represents the boundary between interactions in- 
volving hadronic and partonic degrees of freedom. 
Dotted-lines represent the results of ]5^/. The 
ground state of nuclei is shown as half circle. 



The measured mid-rapidity particle ratios fl,|5T| were fitted with a statistical model [52 
and the results are shown in Fig. 7. The resulting thermal parameters for temperature 
and baryon potential are T c h = 190 ± 20 MeV and [i B = 45 ± 15 MeV, respectively. It is 



not surprising that the temperature parameter is close to the prediction |53|] for RHIC. 

The systematics for chemical parameters is shown in the phase plot in Fig. 8 where the 
dashed lines represent the boundary between interactions involving hadronic and partonic 
degrees of freedom. The dotted lines represent the results of J50| . Results from heavy ion 



collisions at SIS, AGS, and SPS energies are shown as circles. The result from collisions 
at the RHIC energy is shown as star. The lattice prediction on T c (^ 160 MeV) is shown 
at \ib = and the baryon density for neutron star is indicated at the T c h = 0. 
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There are technical caveats in above approaches. The application of a statistical model 
requires that the measurement is done over the whole phase-space, i.e., with 4n particle 
yields. This is because conservation laws that apply to the collisions are only valid for 
global measurements, not locally. In all experiments the coverage of the phase-space is 
limited. Therefore, extracting the yields measured within limited phase-space to 4-7T yields 
is almost impossible, especially for the collider experiments. In addition, we should point 
out that the non-local effect is common to all global conservation laws (baryon number, 
charge, and so on), not only to the strangeness. 

5. Summary 

In summary, the most interesting results from collisions at RHIC are that the system is 
indeed approaching net-baryon free and the transverse expansion is much stronger than 
that from collisions at AGS/SPS energies. 

In order to understand the trend of the collective velocity, an energy scan between y/s NN 
= 20 - 200 GeV, is important. In addition, systematic studies on the anisotropy parameter 
t> 2 and the transverse momentum distributions of 0, Q, and J /ip are necessary as they will 
help in determining whether the collectivity is developed at the partonic stage. 

We are grateful for many enlightening discussions with Drs. P. Braun-Munzinger, 
W. Busza, M. Gyulassy, U. Heinz, B. Holzman, P. Huovinen, B. Jacak, F. Navarra, S. 
Panitkin, K. Redlich, H.G. Ritter, K. Schweda, E.V. Shuryak, R. Snellings, D. Teaney, 
F.Q. Wang, and X.N. Wang. This work has been supported by the U.S. Department of 
Energy under Contract No. DE-AC03-76SF00098. 
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